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ABSTRACT 


The reduction of structural vibrations in ships and submarines is a long standing 
eoncem of the Navy. Wavectide absorbers are very effective devices whicli can be ap- 
plied to this problem. This studv evaluates the increase in vibration damping of a plate 
structure across a broad frequency range using light weight beam waveguide absorbers. 
Viscoelastic and constrained laver beam waveguide absorbers were studied beth theore- 
tically and experimentally. Impedances of the waveguide absorbers at the attachment 
point were predicted using both Bernoulli-Euler and Timoshenko beam theorv for the 
viscoelastic beain and using sixth order beam theory for the constrained laver beam. The 
theoretically predicted impedances were compared with the experimental measurements. 
Results from random vibration tests of a plate structure showed significant increases in 
damping over a broad frequency range (100 Hz - 2 KHz) when the waveguide absorbers 


Were attached on the plate. 
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I. INFRODUCTION 


Suppressing noise and vibrations of ships and submarines is very tmportant tor the 
Navy and has been one of the Navv’s long-standing concerns. [uch work has been done 
on the development of vibration control means such as isolation, detuning. viscoelastic 
damping and dynamic absorption. All of these approaches have been studted rather 
thoroughly including their range of applicability and definite limits on their vibration 
and noise reduction capability within acceptable weight and volume increases. 

Recently. a waveguide absorber concept has been developed [Ref. 1], which can 
provide simple and highly effective vibration control over a wide frequency range. A 
“Waveguide” is a structure along which vibrational waves can travel. If one end of a 
waveguide 1s attached to a vibrating structure, some vibration energy will travel along 
the waveguide in wave forms. If the waveguide is treated with a high energy dissipation 
scheme, the damping of the waveguide causes the ainphtude of the waves to decrease 
as they travel and the waveguide may be expected to remove vibrational energy from the 
structure. 

In previous studies bv Ungar and Kurzweil (Ref 2] and bv Ungar and Wilhiams [Ref 
3], semu-iufimte beams and exponentially tapered semi-infinite beam waveguide absor- 
bers were studied. Tlowever, theoretical prediction of the driving point tmpedances 
showed wide discrepancy from the experimental results since the expertmnents were per- 
formed with finite length beams and the effects of viscoelastic materials on beams were 
not considered in theoretical calewations. 

In this investigation, two kinds of high damping beam waveguide absorber are 
subjected to studv: viscoelastic beam and constrained laver beain waveguide absorbers. 
The driving point impedances are theoretically predicted for the finite length beams of 
lightly viscoelastic behavior. Experimental studies are also performed on the impedances 
of these waveguide specimens and on their contribution to damping increase of a test 


plate structure. 


Wl. THEORETICAL ANALYSIS 


A. LOSS FACTOR CONTRIBUTION OF A WAVEGUIDE ABSORBER 


Consider a linear structure, S, under a harmonic load 
F =/fp sin wt ee 


at location c (Figure | on page 3). The structure $ will vibrate with the enerov sie 
bration 


2 | . 
eae pale i 
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Where Jf represents the total mass of the structure and V,, 1s the magnitude of average 
velocity. The energy dissipated in S per cycle, D, depends on the original loss factor, 4. 


of the structure and the following relation holds: 


Do 
is 2aHh 4 


For a lrghtly damped structure, 7, 1s much smaller than | and D, is much smaller than 
2nit,. Ifa pomt @ of a second body B. a waveguide absorber, 1s attached to body 5 
at Location b of body S . internal force, F,, will interact between body S and body #@ and 
some vibration energy, D will be dissrpated mn body B (Figure I on page 3). Moreover, 
this internal force influences the vibratronal motion of body S and mav reduce the en- 
ergy of vibration. HV, of body S under the same harmonic load at the same location ¢. 
The interaction between bodies S and B depends on the impedances of bodies S and 
Bat the attachment point, 2, and Z,. Impedance is the ratio of the harmonic force, F,. 


acting Oa pomt @ of the structure to the velocity, J%, of the point. 


I 
Ze {2.4} 
J 
and can be deserrbd bv a complex number, 
Z=R+iN (2) 


to 





Figure 1. Vibrating structure with waveguide absorber. 


where R is the real part of Z and resistance component, or simply resistance, and .¥ is 
the imaginary part of Z and reactive component. or reactance. 
The energy loss per cycle, D, of a structure by an attached waveguide absorber ts, 


as shown by the previous study [Ref. 2. and Ref. 3], 
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where w denotes the vibration frequency in radian'sec. V, represents the velocity of the 
structure at the attachment point of the structure before the attachment of the wave- 
guide absorbers, Z, and Z, denote the impedance of the waveguide absorber and the 
structure at the attachment point, and 2, represents the real part of the impedance Z,. 
From the equations 2.2, 2.3 and 2.6, the energy loss factor. #7, due to the 


contribution of an attached waveguide absorber can be expressed as follows: 
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Equation (2.6) shows that the vibration energy loss of a structure through the wav- 
eguide absorber depends on the magnitude of the velocity, 1, at the attachment point 
on the structure. If this attachment point is a node for which }, = 0. then the energy loss 


factor. 7. equals zero. If the attaelinent point is an antinode, then | > can be rela= 
tivelv large. Equation (2.7) indicates that the loss factor contribution is small if the im- 
pedance of the beam, Z, , 1s considerably greater than the mmpedance of the strircttme: 
Z,.. If Z, is much smaller than Z,, then the structure will not be alfected by the absoroee 
and the vibration energy of the structure, I, in equation (2.2), wili not be redueed mueh. 
Therefore, iunpedance matching between the structure and the waveguide absorber ts 
uportant to get the highest energy loss factor, 7. So. analytic methods whieh ean predict 
the impedance of beams using the theorv of elastic Wave propagation in beams are de- 


veloped in the next section. 


B. ISEPEDANCE OF BEAM WAVEGUIDE ABSORBERS FROM WAVE 
PROPAGATION THEORY 
I. VISCOELASTIC BEAN USING BERNOULLI-EULER BEAM THEORY 

Previous research [Ref. 3] considers the impedances of the infinite viscoelastic 
beam using Bernoulli-Euler beam theorv. In this section, this previous work is expanded 
to include the finite viscoelastic beam. The viscoelastic beam will be exerted harmonically 
at its center bv a transverse force, F as shown tn Figure 2 on page 6. For a beam which 
is free of lateral loading and under an assumption that cross-sectional areas remain plane 


and normal to the neutral axis, the equation of motion becames [Ref. 2]. 
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where (x7) 1s the transverse displacement, and E* is the complex moduli Gf the beam 
E* = E(1 + imp). (2.9) 


and, 


i: modulus of elasticity of the beam 


Te lOSS tecuormo. (he VISCOClMilemmarcrial 

m : mass per unit length of the beam 

Deedes more Cea maternai 

Fl MChOsceccemotial area of the beat 

[: moment of inertia of the cross sectional area 

The general form of +(x.) for a propagating harmonic wave in this Bernoulli- 


Euler beam ts 
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Where 4* 1s the complex wave number, and q is the frequency of the propagating wave. 
Inserting the wave equation (2.10) into the equation of motion (equation 2.8), the dis- 


persion relation 1s obtained as follow: 
E*(k*) = pAw* (2.11) 


This equation gives the four different complex wave numbers as function of frequencies, 


| 1 
x : Sz tan , nf tails 4. 
AX = ik(/1 +92) Meos( =) — i sin( =) (2.12.4) 
ae tan y tan 7 
x : / De veaeat c ve > ae ve 
ky = —ik(./ 1 +2)" [cos( ae el rae (2.12.5) 
= tan” |y tan 17 
ky =— K(f} + ns, 7 Feos( > — / sin( aa ) (Ae) 
a tan y tan 7 
k= kf 1 + 5,7 Leos ae — i sin{ er ) (2.12.<) 


where subscript ve 1s an abbreviation of Viscoelastic beam using Bernoulli-Euler beam 


theorv and & 1s the wave number, given by 
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meeemotes the radius of gvration of the cross sectional area, ( V/ aed and c, denotes the 
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longitudinal wavespeed im the beam material (om jpineretone, wie ceuniion (2-14) 1s 


expanded using wave number ky, kx, Ay and Ax. 





Figure 2. Coordinates and sign conventions for a semi-infinite (finite) viscoelastic 


beam. 


Ge) = (eer x zi ieee 3 ¥,e"% AS wie Koel ema (2.1) 
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Where Y,, ¥;, ¥; and ¥, are constants whose values depend on the boundary conditions 
and x is the axial coordinate. The complex transverse force, F, is denoted by 
dt s{x) 


a 
ciX 


F(x) = E*I (2.15) 
a. INFINITE VISCOELASTI@ 3 iy 

Let’s consider the infinite viscoelastic beam. If yx) is to remain finite for a 

large x, then ¥, must vanish. If there is no wave comung toward the origin from the po- 

sitive x direction, then ¥, must vanish. Y; and Yj; in equation (2.14) are determined from 


two boundarv conditions at x = 0, as functions of vy and @, as follows: 


y(O) =y (2.1Gie 


Av{O) = 
———- = @ (25) 
ax 
At the driving point, the input transverse force, /(0), can be calculated using 
the equation (2.15), which can be expressed using impedance and velocity at « =0 as 


follows: 
F(O) = Z, V(O) + Z,Q(0) (2.18) 


The complex transverse velocity, V, and the complex angular velocity, ©. are related to 


the corresponding displacement. y and @ of equations (2.16 and 2.17). as follows: 
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Z, represents the tmpedances of the transverse force at the driving point of the center 
of the infinite viscoelastic beam. If there 1s no rotation at the center of the beam. that 


ieee.) = 0, tlle equation (2.18) becomes 
F(O) = Z, 1'(0) 2.20) 


Therefore, the impedance of an infinite Bernoulli-Euler beam at the driving point. Z.,,.. 


can be calculated from 


FiO) 1 


iQ) J 


where subscript ive is an abbreviation of the Infinite Viscoelastic beam using Bernoul- 
li-Euler beam theorv. 
Demet VASCOELASTIC BEAM 
Por auiniielopeam, (Our constants ),. 14, 1%, and ¥,, are determined (rom four 


weundary conditions at x = 0 and x =/, as follows: 


3(O) = ¥ ine 
a = (2.23) 
QX 
Py 
See #0 = () (2.2.4) 
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From the same considerations as those in the previous seetion, the impe- 
dance Z,,, 
Fr) 


ey 
Zyl) ay io a (26) 





1s determined as function of frequencies. Where the subscript fre 1s the abbreviation of 
the Finite Viscoelastic beam using Bernoulli-Euler beam theorv. 

2. FINITE VISCOELASTIC BEAM USING TIMOSHENKO BEAN THEORY 

The Bernoulli-Euler beam theorv was obtained using two assumptions which 

neglect shear deformation and rotational inertia. The Bernoulli-Euler beam theory aia 

be useful for low frequencies and long wavelengths. However, the Timoshenko beam 

theory gives more accurate simulation at wider frequency ranges and wave lensthis: 

Therefore we will analvze the impedances of the finite viscoelastic beam using Tinosh- 

enko beam theory in this section. Two tndependent variables, the transverse displace- 

nent. 3(x.2), and angular displacement, @(x.2), are defined in Timoshenko beam theorv. 


Beaim motions are governed bv the following equations [Ref. 4]. 
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Where subscript rt is an abbreviation of Viscoelastic beum using Timoshenko beam the- 
ory, and 

G ; shear modulus of the beain 

E: elastic modulus of the beam 


#,.; loss factor of the viscoelastic material 


a; cross sectional area of the beam 

I: moment of the inertia of the cross sectional area 

po: densitv of the beam material 

«x : Shear deformation factor 

The general forms y(x.z) and O(x,2) for the propagating harmonic wave in a [1- 


moshenko beam are 


eee oN) 
y 0 
b(x.0) = Dee no! (2.32) 


where k< is the complex wave number, and w 1s the frequency of the propagating wave. 


Inserting the wave equations (2.31 and 2.32) into the governing beam eyuations (2.27 
and 2.28) becomes 
, p iq : 2 r ° a ,4 
(wo* — as(k*)*)¥y + ( —iafk™) Dy = 0 (2.33) 
- I a2 x \2 Pe 
(achele W873 om = aalke STA yy Se (2.341) 
The disperston relation of waves 1s obtained from equations (2.33 and 2.34) 
i) ae DG alan) Ao ee: 
avas(k*) —(ay + as)w(k*) +(@ — apkjw’) =0 (2.35) 
with the relationship between 1, and ®),. 
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The equation (2.35) gives the four different wave numbers as function of wave frequen- 


cies, 
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Therefore, the equations (2.31 and 2.32) are expanded wsiie A*ha5 sce 
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The equation (2.40) becomes the following equation 
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(2.42) 


For a finite viscoelastic beam with a length /, which is free at x =/ and excited 


at x = 0, four boundary conditions are defined as: 
JO) =y 
PO) = ¢ 
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(2.-16) 


froin those boundarv conditions, the constants 1, ¥,. 1, and ¥, are deternuned 


as functions of y and @. At the driving point, x=Q0, input 
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transverse {ore 


F(0) = — S(O), can be calculated using equation (2.46), which can be expressed using 


impedances and velocity at the driving point as follows: 
F(O) = Z, F(O) + Z,0(0) (2.47) 


where the complex transverse velocity, ’, and the complex angular velocity. <D, are re- 


lated to the corresponding displacement as follows: 
5 (2.48) 
If there is no angular displacement, @, at the center of the beam, that 1s, D(O) =, the 
equation (2.47) becomes 
FO) = Z, ¥(0) (2,19) 


Therefore, the unpedance of the Timoshenko beam at the driving point, Z,,, can be cal- 


culated using the following equation. 
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where subscript fvr is an abbreviation of the Finite Viscoelastic beam using Timoshenko 
beam theorv. 
3. FINITE CONSTRAINED LAYER BEAMNI 

In this section, impedances of the constrained laver beam, which has viscoelastic 
material between the two elastic lavers, are studied. The theory of a damped sandwich 
beam was developed bv R.A. Ditaranto [Ref. 5], who extended Kerwin’s [Ref. 6 | bas- 
ically similar analvsis. D.J. Mead and S. Markus [Ref. 7] expanded the concept bv con- 
sidering these earlier works. The result of their efforts 1s a sixth order differential 
equation of motion which ts expressed in terms of the transverse displacement, 1, and 
longitudinal displacement, u. for the constrained liver beam. The constrained laver 
beam will be excited harmonically at its center bv a shear force, S, as shown in 
Figure 3 on page 12. This equation of motion leads to the impedances of the sandwich 
beams with six boundry conditions, which are used for evaluating impedances depending 
on frequency, viscoelastic material properity, and two elastic material properties. 


From the [Ref. 7], the governing sixth order differential equation is 
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Figure 3. Coordinates and sign conventions for a semi-finite constrained layer beam. 
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with the relationship between +(.v.2) and u(x.) expressed as: 
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The parameters in this equation are defined as follows: 





Cae 9h 


Sa a ares 





Gay FE H,E,H ; 

Z = — (| (2.34.5) 
Da ao eee TT 

D,= BoE aE E31, (2eSsec) 
C,= Hy, +—(H, +H) (2.54.cf) 
Hy . 
Le 3 (2.54.e) 
Hy; ie, 
I, ae (a) 
G*=G(1+ i, .) 


where subscript c is an abbreviation of constrained laver beam, and 
G : shear modulus of the viscoelastic laver 
E : elastic modulus of the elastic layer 
jaooss wactor of the viscoelastic laver 
{: moment of the inertia of the cross sectional area 
ff ; thickness of the constrained laver beam 
nt: mass per unit length of the whole three layer section. 
The general form of transverse displacement, 3(.x.2), and longitudinal displace- 


ment, u(x,z), for the propagating harmonic waves in a constrained layer beam are 


i(k” x—at) 


(xt) = Toe (50) 
con ten tr (2.57) 


where &* is a complex wave number and w 1s the frequency of the propagating wave. 
After inserting the wave equations (2.56 and 2.57) into the sixth order beaim equation 


(2.51), the dispersion relation is obtained as follows: 
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There are six different roots with six different complex wave nuimberssimetiemaees 
equation, which means a wave of six different modes propagate through the constrained 


laver beam. From the equations (2.56 and 2.57), we iaytGairance Calera eens 
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The general form of the transverse displacement equation (2.56), v. and longi- 
tudinal displacement equation (2.57). wu, are obtained using the six complex wave num. 


bers and a ratio K,. The equations (2.56, 2757 and) 259 eas ae 
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n= 1, 2,3,.4,5. 6 For a constrained laver beam with a finite length / which isis 
at x =/ and excited at one end at x =0, the six boundary conditions are deimiedsasmiae 
1OAvS |Inera |: 
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litieesimeconstants }. (7 =i) 2) $e4e sand 6) are deterinined (rom tliese six boun- 
darv conditions and are a function of y, 0 and uw. At the driving point, x =9, the input 
transverse force. /(0) = — S(0). moment. 3/(0)}, and normal force, \(0), can be calculated 
using equations (2.67, 2.68 and 2.69). These input force and moment can be expressed 


using impedances and velocities and angular velocity at the driving point as follows: 
F(O) = Z, VQ) + Z,Q(0) + Z,U(0) (80) 


Where the complex transverse velocity, }’, and the complex angular velocity, 
Q. and the complex longitudinal velocity, U’, are related to the corresponding displace- 


ment as follows: 
V ; on o 
ee ey ee et) 


If there is no rotational displacement, @. and no longitudinal displacement. uv. at the 
Sennen of the beam, Z, represents the impedances of the transverse force at the driving 
fma@eeo! the constraimed layer beam. That is, €(0)=0 and U(G)=0. Therefore. the 


equation (2.70) becomes 
F(O) = Z, ¥(0) 2) 


Therefore, the impedances of the constrained layer beam at the driving point, Z,, can 


be calculated using the following equation. 


FAO) ila 


—1ny S$ 


Z {O) = 


Where subscript c is an abbreviation of the constrained laver beam. 


C. COMPUTER IMPLEMENTATION 
The computer programs using Reduce and FORTRAN languages are used for cal- 
culation of the theoretical impedances from the beam theories in this paper. Reduce was 


developed by Anthony C. Hearn [Ref. 8], and is a programming language which solves 


I 


algebraic operations non-numerically. It can manipulate polynomials in a varity of forms 
(substitution, differentiation and integration) and solves une or more sirnultancous al- 
gebraic equations. After the impedances were evaluated by a Reduce program with 
constant variables. a FORTRAN program calculated the impedance with real dimen- 
sions and properties of the viscoelastic and constrained laver beam waveguide absorber 
for the 0-2000HZ frequency range. All constant variables of the Reduce program were 
defined in each FORTRAN program. We tried to use the same notations for the Reduce 
program. the FORTRAN program and the text (Appendix A). 

Three Reduce programs and a FORTRAN program were used to study cach theory. 
Thre first Reduce program defined the moment, force and displacements which icluded 
the constants }¥,., where # equals | and 3 for an infinite viscoelastic beam; 7 equals |, 2, 
3 and 7 for a finite viscoelastic beam and » equals 1, 2, 3, 4, 5, and 6 for a constraimed 
layer beam. ¥, can be solved with boundary conditions. After running the first Reduce 
program, neW constant terms are defined for ¥,. The second Reduce prograin includes 
those new constant vanables which were the results of the first Reduce program. [he 
second Reduce program gives the constant variables Y,. Finally, the third Reduce pro- 
grain can evaluate the impedance of the beam with longitudinal and angular displace- 
ment (and normal displacement for the constrained layer beam). The process of the 
constramed layer beam is similar up to the second Reduce program. However, the 
computer finds 1t hard to handle 6 simultaneous equations with 6 non-numerical vari- 
ables, ¥. Therefore. we solved ¥,, ¥. and }, in the second Reduce program) andegaew 
substitute these values to ¥,, }; and ¥, in the third Reduce program. 

Wave numbers as functions of frequencies, i.e. dispersion relation, of waves in way- 
eguide absorbers were described in previous sections. It was easy to find the wave num- 
ber of a viscoelastic beam through Bernoulli-Euler beam theory and Timoshenko beam 
theory in equations (2.12 and 2.38). However. wave numbers of the constrained laver 
beam were expressed with a complex six order polynomial. Therefore, they were calcu- 
lated using Newton's method. The equation (2.58) 1s rearranged using real constants 
bh 


c, and d, and mraginary constants 4, and d. 


r* 


4 
a 


(hay +(b, + ib) RK) + c(K") + (d, + id) = 0 (2.74) 


where 


> 
x 
i 
—~ 
on 
x 
~~ 
i) 
~~ 
i) 
~~) 
tan 
S 
~ee” 
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hb, = [17 Dy + BZ) (pe) G) 
pa =o” (Pema) 
; D ead 
a 
/ — By (2 as ) 
d. = Re( —————_ aise 
r D, 
an SSiantey: (751) 
a= 11 : fen 
i oS 


Let the complex variable K* = 1+ ig (t and q are real numbers). The equation (2.74) 


becomes two equations coupled with two variables (¢ and q). 
emia emo bh ( —3¢ — 2bg Pere bg + d,) (2.76) 
etigv=@ + bq? + (30° — 2b,t — eq + (—bt° — d) (2.77) 


Starting from initial guesses (f,, g,) Which are near the roots iteration was continued until 
either the successive ¢ and g values converged to the certain values or the value of the 
function is sufficiently near zero. After a few iterations of equations (2.76 and 2.77) 
using a FORTRAN program (Appendix A), we can get three different square of wave 
number values. A*, were obtained. So, finally. this method calculates six diflerent wave 
meenioers, kx, from + JkKé, mules # Clas ih 2.1)... 5, and 6. 

Impedances of a 20” viscoelastic beam were calculated using these computer pro- 
grams for viscoelastic beam properties of 200,000 psi constant shear modulus and 0.2 
energy loss factor. The results of these sumulations are shown in [igure + on page 19 
and Figure 5 on page 20 which reflect sharp resonance frequencies. These sharp reso- 
nance phenomena reduced for the same beam with higher energy loss factor of 0.5. For 
Gases Using the Bernoulli-Euler beam theory, the impedances of the finite viscoelastic 
beam approached to the infinite viscoelastic beam waveguide absorber with increasing 
energv loss factor and with increasing frequency. The impedance using Timoshenko 
beam theory differs from the impedance from the Bernoulli-Euler beam theory 
(Figure 6 on page 21 and Figure 7 on page 22), because Timoshenko beam theory 


considers shear deformation and rotary inertia effect of the viscoelastic beain. 
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Impedances of a 20” constrained laver beam with viscoelastic layer properties of 200 psi 
constant shear modulus and 0.5 and 1.0 energy loss factors were calculated. Figure 8 
on page 23 and Figure 9 on page 24 compare the impedances vs. frequency depending 
on the energy loss factor of the constrained layer beam. The tmpecance of the high 
dumping constrained laver beam shows the smooth and shifted resonance frequencies 


compared with the low damping constrained laver beam. 
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Figure 4. Real part of the 16” viscoelastic beam waveguide absorber theoretical im- 


pedances (7 = 0.2) . 


19 


1400 1800 £800 2000 


1200 





600 6800 1000 
FREQUENCY (HZ) 


IMPEDANCES VS. FREQUENCY 


400 


ST o'r s‘0 0°0 s'0- 
(NI/D9S—4AG1) ‘GAd WI AUVNIDVAI 


oe 0°2 


Figure 3. Imaginary part of the 16” viscoelastic beam waveguide absorber theore- 
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tical tmpedances (17 = 0.2) . 
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Figure 6. Real part of the 16” viscoelastic beam waveguide absorber theoretical im- 


pedances (7 = 0.5) . 
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Figure 7. Imaginary part of the 16” viscoelastic beam waveguide absorber theore- 


tical impedances (7 = 0.5) . 
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Figure 8. Real part of the 20” constrained layer beam waveguide absorber theore- 


tical impedances. 
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Figure 9. imaginary part of the 20” constrained laver beam waveguide absorber 


theoretical impedances. 


DEX BRENT 


A. IMPEDANCES OF THE TEST PLATE 

A rectangular aluminurn plate with clamped boundary condition was sclected as the 
test structure, because it has enough closelv space modes tn the frequency range 
(100Hz-2Kz). An aluminum block (30” x 22” x 2”) was carved inside to make a 24" x 16° 
x 0.3125" plate with clamp boundary condition. 

The impedances of the test plate were measured from impact hammer tests. The 
aluminum plate was excited by a PCB 086B03 impact hanimer. Input impact force was 
measured bv a force transducer at the tip of the impact hammer. Responscs were meas- 
ured bv an Endevco 2250A-10 piezoelectric accelerometer which was attached under the 
linpact point for locations 1, 2 and 3 (Figure 10 on page 26 and Figure Il on page 
27). Location | is the center of the plate. Location 2 is the symmetric point of location 
4, where a shaker will be attached during the damping measurement tests. Location 3 ts 
a general point which is not on anv line of symmetry. A baseband measurement for each 
location was taken from 0 to 2000Hz. An average of 15 data samples was taken for each 
measurement and analvzed using a HP-3562A dvnamuc signal analyzer. Impedance of 
the aluminum plate at each location was determined atter a predetermined data 


processing procedure using the HP-3562A dynamic signal analyzer. 
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Figure 10. Arrangement for test plate impedance measurement. 
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Figure lJ. Test plate impedance measurement configuration. 


B. IM[PEDANCES OF THE WAVEGUIDE ABSORDEE 

Two viscoelastic beams and two constrained laver beams were selected as the wav- 
eguide absorber specimens. Dimensions and phvsical properties are shown in Table 1. 
Mechanical properties of a viscoelastic material depend on temperature and vibration 
frequencies. The loss factors and shear modulli of the viscoelastic material used in the 
specimens at room temperature (75°F) are shown m Figure 12 on pace 50a 
Figure 13 on page 31. These data were obtained from Design Data Sheets in a book by 
Nashif, et. al. [Ref. 9]. Analytic expressions of the loss factors and shear modult were 
determined from curve fits to these data. The following four equations represent the loss 


factors and shear moduli of the LD-400 and ISD-112 viscoelastic material. 


ry = 0,65 ¢170527321 tow Fy 1 3.5 
aieo eo! log( =~) |! 3 Gal 

G, = 0.00002503/" —0.1752f" +457.5883f +29280 (3.3) 
G, = —0,00001258¢° +0,2472f + 74.988 (3.4) 


f represents the frequency. 7, and y, represent the loss factor of LD--100 and ISD-1}2 
viscoelastic material. G, and G, represent the shear inodulus of LD-400 and [SD-1H12 
viscoelastic material. Figure 14 on page 32 shows waveguide absorbers which were used 
in this study. 

[mpedances of waveguide absorbers at the attachment points were determined from 
steady state random vibration tests. The waveguide absorber was excited by a Wilcoxon 
F4 F7 vibration generator using a band limited white noise signal (O-2000Hz) from a 
HP-3562A dynainic signal analyzer (Figure 15 on page 33). The waveguide absorber was 
mounted on the Wilcoxon F4. 7 vibration generator (Figure 16 on page 34). Force and 
acceleration were measured using a force transducer and a piezoelectric accelerometer in 
the impedance head of the Wilcoxon F4.F7 shaker. These signals were analyzed 
HP-3562A dynamic signal analyzer. 10U0 data samples were averaged for each measure- 
ment. 

The measurement of rmpedances included the impedances of bolts and nuts and 
connecting alununum pieces. Therefore, the impedances of bolts, nuts, and connecting 


aluminum preces, the so-called mass effect, were measured separately. Total impedances 


minus mass effect represented the impedances of the waveguide absorbers. These 


experimental results are compared with the theoretical predictions in Chapter IY. 


Table 1. DIMENStONS AND PROPERTIES OF THE WAVEGUIDE ABSOR- 
BER SPECIMENS. 
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Figure 14. Waveguide absorbers. 
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Figure 15. Arrangement for waveguide absorber impedance measurement. 
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Figure 16. 


Wilcoxon F4/F7 shaker with waveguide absorber mounted. 


C. DAMPING OF THE PLATE 

Modal damping values of the plate were determined from structural transfer func- 
tions at a driving point and between input and output points on the test piate which 
were measured using randoin excitations. The HP-3562A dynamic signal analyzer was 
used in the measurement of the transfer functions. The Wilcoxon F3 vibration generator 
excited the plate at location 4 using random oscillation signals from a HP-3562.A dy- 
namic signal analvzer. Input forces and output responses at the driving point were 
measured using a force transducer and a piezoelectric accelerometer in the impedances 
head of Wilcoxon F3 shaker (Figure 17 on page 36). The damping measurement con- 
figuration is shown in Figure 18 on page 37. It was necessary to measure responses at 
several points other than the driving point to accurately measure damping values of 
closelv spaced modes. These were measured using the Endevco 2250A-10 piezoclectric 
accelerometer. Modal damping values of all the modes in the frequency range of 
100-2000Hz were deterinined from the response functions using the Half power Band 
Width method [Ref 10] and a curvefit method (Ref. 11] 

Modal damping values of the plate without a waveguide absorber were measured 
as the baseline. A waveguide absorber was attached to the plate at location 1. [he 
waveguide absorber was changed after each measurement of the damping for four dif- 
ferent waveguide absorbers which were described in section 4-B. After the damping was 
measured at location 1, measurements were taken at location 2 and location 3. 

The most effective waveguide absorber at each location was selected. The 20° vis- 
coelastic beam waveguide absorber, the 16” constrained laver beam waveguide absorber 
and the 16” viscoelastic beam waveguide absorber were most effective at location 1, 2 
and 3, respectively. Therefore, the damping with two waveguide absorbers attached at 
locations | and 2, and the damping with three waveguide absorbers at locations 1. 2 and 


3 were also measure, to see the effect of multiple waveguide absorbers. 
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Figure 17. Arrangement for test plate damping measurement. 
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Figure 18. Damping measurement configuration. 
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[V. RESULTS AND DISG@USsIOrs 


A. INIPEDANCES OR Tithe ee lia 

Impedances of the test plate were meusured from input hammer tests as described 
in Section 3.A. Since the aluminum test plate is a hghtly damped structure, the impe- 
dances at resonance frequencies are very low and at anti-resonance [frequencies are verv 
high. Phases with 1809 phase shilt are near + 909 both at the resonance tide 
antui-resonance frequencies. 

Since location | is the center point of the plate. only small numbers of resonance 
and anti-resonance [requencies appear (Figure 19 on page 39 and Figure 20 on page 
40). Neither location 2 nor location 3 are on any line of symmetry of the test ¢laie: all 
the resonance and anti-resonance points appear in the mnpedance vs. frequency dia- 
grams. (Figure 21 on page 41, Figure 2? on page 42. Figure 23 On pace eis 
Figure 24 on page 44). However, the anti-resonance frequencies at location 2 and those 


at location 3 are different. 


Figure 19. 
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Impedance magnitude of test plate at location 1. 


Figure 20. 
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Figure 21. 
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Impedance magnitude of test plate at location 2. 


4] 





LOCATION2 Dk 


H=z 


Fxd Y O 


‘oure 22 
htouresc2: 


Ex/Fre 


a 
> 

O 

x 

Q 

a 

> 

G 

ts 

0. 

) 

lJ 

X 

CG 

Lu 

XO eae 0 

iL © os Q) 

ro N Oo 0 0) 
2 Se C g 
ab 1 QO 


Impedance phase of test plate at location 2. 








oe 1) 





LOCATION2 


| es 


Bo Glan 


Figure 23. 






LOCATIONS 


fl 

U) N 
W 1G 
Y 

Q O 
LJ 

YO a > Bos 
ly 2 ; “al 

ne) 2s! Oo 
0) ih wen a MX 
ale U eau 


Impedance magnitude of test plate at location 3. 


43 


Figure 24. 
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B. IMPEDANCES OF THE WAVEGUIDE ABSORBER 

Two different kinds of waveguide absorbers, a viscoelastic beam and a constrained 
laver beam. are investigated in this study. In a previous study by Ungar and Wilhams 
[Ref. 3]. tmpedances were predicted only for the infinite viscoelastic beam using Ber- 
noulli-Euler beam theory and the theoretical prediction showed wide difference from the 
experimental results. In this study. this previous work is extended to include the impe- 
dances for the finite viscoelastic beam. Since Bernoulli-Euler beam theory does not in- 
clude the shear deformation and rotary inertia effect of the viscoelastic beam, the 
predictions from tle Timoshenko beam theory are studied to investigate these effects on 
Wave propagation and impedances, and are compared with the results from the Ber- 
noulli-Euler beam theorv. The impedances of the constrained laver beam waveguide ab- 
sorber were predicted using the sixth order beam theory. 

The impedances from these theoretical predictions are compared with the impe- 
dances from the experimental measurements for the 16” and 20” length waveguide ab- 
sorbers of two different kinds. In the calculation of theoretical predictions, ithe loss 
factors and shear moduli of viscoelastic material were varied as function of frequency 
as described in Section 3.B. During the experiment. the Wilcoxon F4. F7 vibration gen- 
erator produced characteristic errors which occured at different frequencies depending 
on the mass of the waveguide absorbers. Therefore the viscoelastic beam has unreliable 
data around 70QHzZ and the constrained laver beam has unrehable data around 350FIz. 

Figure 25 on page 47 and Figure 26 on page 48 show the three theoretical impe- 
dances compared with experimental data. The Bernoulh-Euler beain theory for finite 
beams shows some standing wave effect at low frequency range due to reflection at the 
free end of the beam, but approaches the results of the infinite beam as frequency in- 
creases. The results of the Bernoulli-Euler beam theory do not quite correspond to ex- 
perunental results for the viscoelastic beam. Iimpedances from the Timoshenko beam 
theory closely follow the experimental results for viscoelastic beams as shown | in 
Figure 25 on page 47 and Figure 26 on page 48. These show that the shear deformation 
and rotarv inertia effects are significant for the impedance prediction of beam type 
waveguide absorbers. 

Figure 27 on page 49 and Figure 28 on page 50 represent the real and imaginary 
parts of the 16” constrained laver beam waveguide absorber and Figure 29 on page S| 
and Figure 30 on page 52 represent the real and imaginary parts of the 20” constrained 


layer beam waveguide absorber. In both beams, the theoretical results of the impedances 


45 


using the sixth order beam theorv and the resonance frequency patterns are almost the 
same as the experimental result. 

Both Timoshenko beam theorv and the sixth order beam theory predict resonance 
frequencies shifted amount from the experimental results as frequency increases, and the 
Impedance magnitudes of experimental results are bigger than the theoretical results. 
These difference between theoretical and experimental impedances may derive from the 
inaccuracy of shear modulus and loss factors which values were selected {rom 73°F room 
temperature and changed depending on frequency. Hlowever, the theoretical prediction 
using Timoshenko beam theorv for the viscoelastic beam and the sixth order beam the- 
ory for the constramed laver beam follow the same trend (figure 12 on page 30 and 
Figure 13 on page 31). 

[impedances of waveguide absorbers are very low and do not show anv sharp vari- 
ation for all frequencies in the range tested. which are quite different from the inpedance 


of the test plate of very lightly damped structures. 
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Figure 25. Real part of the 16” viscoelastic beam waveguide absorber impedances 


at the center of the beam. 
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Figure 26. [maginary part of the 16” viscoelastic beam waveguide absorber iimpe- 


dances at the center of beatin. 
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Figure 27. Real part of the 16” constrained layer beam waveguide absorber impe- 


dances at the center of the beam. 
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Figure 28. Imaginary part of the 16” constrained layer beam waveguide absorber 


impedances at the center of bean. 
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Figure 29. Real part of the 20” constrained layer beam waveguide absorber impe- 


dances at the center of the beam. 
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Figure 30. Imaginary part of the 20” constrained !aver beam waveguide absorber 


impedances at the center of beam. 
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oe DAMPING OF THE PLATE 

The purpose of this experiment is to investigate the effect of light weight waveguide 
absorbers on damping increase tin plate type structures over a wide frequency range. 
Damping of the test plate structure was measured with an attachment of each one of 
four different waveguide absorbers (16” or 20” viscoelastic or constrained laver beam) 
at three different locations. The frequency response and the moda] damping values of the 
aluminum test plate with a waveguide absorber were compared with those of the same 
plate without anv waveguide absorber over a wide frequency range (0-2000HZ). 

The modal damping values of the test plate without a waveguide absorber are shown 
in Figure 31 on page 56 and Figure 32 on page 57. They show that more than 2U lightly 
damped modes are m the frequeucy range between 1G0f{[z to 2000HZ with the lowest 
imequency 195EIz. 

The damping contribution of the same Waveguide absorber to a same plate was dif- 
ferent depending on the location of the waveguide absorber attachment poit. For the 
20” viscoelastic beam waveguide absorber, there was marked damping increase for 
modes lower than 1100Hiz when it was attached at location | (Figure 33 on page 58 and 
Figure 34 on page 59). However, it was more effective for modes higher than 1200Hz 
When it was attached at location 2 (Figure 35 on page 60 and Figure 36 on page 61). 
When it was attached at location 3 its contribution to plate damping increase was el- 
fective for wider frequency range of 25011z-1500Hz (Figure 37 on page 62 and 
Figure 38 on page 63). This trend applies to other waveguide absorbers, though the 
magnitude of damping contribution and the frequency ranges were different for cach 
case (Figures in “Appendix B). 

Different modal damping effects with a different waveguide absorber at the same 
location were also observed. At location I, the 20” viscoelastic beam waveguide absorber 
Was inost effective for the frequency range lower than 1lOQHz (Figure 33 on page 38 
and Figure 34 on page 59) and the 16” constrained laver beam waveguide absorber was 
most effective for the frequency range higher than 1200I1z (Figure 39 on page 64 and 
Figure 40 on page 65). At location 2, the 20” and 16” viscoelastic beam waveguide ab- 
sorbers showed moderate contribution to plate damping and the 2U" constrained laver 
beam waveguide absorber contributed the least amount to plate damping. For frequency 
ranges over 1200Hz, the 16” constrained layer bean waveguide absorber produced the 
largest damping increase. At location 3, the 20” and 16” viscoelastic beam waveguide 


absorbers and the 20” constramed layer beam waveguide absorber produced moderate 


o> 


damping ticrease over a wide frequency range of 250Hz-1500HZ. However, the 167 
constrained laver beam waveguide absorber showed the least damping increase. 

After the selection of optimum damping at each location, two wave guide absorbers 
(20” viscoelastic and 16° constrained layer beain) were attached at locations 1 aiidezeam 
three waveguide absorbers (20° viscoelastic, 16° constrained layer and 16 “ViseGelamam 
beam) were attached at locations {, 2 and 3 to see the effect of multiple waveguide ab- 
sorbers on the test plate damping (Figure 41 on page 66). Figure 42 on page 67 and 
Figure 43 on page 68 represent frequency response and modal damping of the plate with 
two waveguide absorbers. Figure 4J4 on page 69 and Figure 45 on page 70 represent 
frequency response and modal damping of the plate with three waveguide absorbers. 
Damping of the plate mncreased with increasing number of the waveguide absorbers 
(Figure 42 on page 67 and Figure 44 on page 69). The damping increase resulting from 
attachment of two or more absorbers resemble the sum of the seperate result of each. 

Impedance of the test plate at location | 1s compared with the driving point impe- 
dance of the 20° viscoelastic beam in Figure 46 on page 71. Thus figure shows that at 
resonance frequencres the impedance of the test plate ts much smaller than that of the 
waveguide absorber. Figure 47 on page /2 and Figure 48 on page 73 show the com- 
parsion between the rmpedances of the test plate at location 2 and location 3 with the 
drving point of the 16” constrained laver beam and 16” viscoelastic beam, respectively. 
These figures show that at low frequency range (<4OOFIzZ) the impedances of the test 
plate were higher than those of waveguide absorbers and at high frequency range 
(>400HZz) impedances of waveguide absorbers were higher at some resonance ftrequen- 
GICs: 

From equation (2.7), the loss factor contribution, 7, of a waveguide absorber ts 


proportional to the following 3 terms: 








R, (4.1) 
! (42) 
= ee 

+I 

fio . 
ma (4.3) 
ny 


The first term is determined from waveguide characteristics only and the second term ts 
determined from the impedance ratio between the test structure and a waveguide ab- 
sorber at the attachment point. However, the third term depends on the interaction be- 
tween the test structure and a waveguide absorber and the whole response of the test 
plate. If Z, is much smaller than Z,, then F, will remain constant before and after the 
attachment of a waveguide absorber. Previous studies [Ref. 2 and Ref. 3] were based 
On this assumption and gave musleading indications that expernnental loss factors would 
be much higher than theoretical predictions. In the present study. as shown in 
Figure 46 on page 71. Figure 47 on page 72 and [igure 48 on page 73, the impedances 
of the test plate were similar in magnitude to impedances of waveguide absorbers at re- 
sonance frequencies and the loss factor contribution due to waveguide absorbers could 


not be predicted using equation (2.7). 
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Figure 34. The driving point frequency response of the test plate without a wave- 


guide absorber. 
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Figure Modal damping factors vs. frequency of the test plate without a wave- 


guide absorber. 
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Figure 33. The driving point frequency response of the test plate with a 20” viscoe- 


lastic beam waveguide absorber (dashed) at location [| and without 


(solid). 
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Figure 34. = Mlodal damping factors vs. frequency of the test plate with a 20”  yis- 


coelastic bean waveguide absorber at location 1 and without. 
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The driving point frequency response of the test plate with a 20” viscoe- 


sy 


lastic beam waveguide absorber (dashed) at location 2 and without 


(solid). 
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Figure 36. Mfodal damping factors vs. frequency of the test plate with a 20” yiscue- 


lastic beain waveguide absorber at location 2 and without. 
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Figure 37. 
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The driving point frequency response of test test plate with a 2 0” vis- 


coelastic beam wayeguide absorber (dashed) at location 3 


(solid). 
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Figure 38. Modal damping factors vs. frequency of the test plate with a 20” viscne- 


lastic beam waveguide absorber at location 3 and without. 
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Figure 39. The driving point frequency respouse of the test plate with a 16” con- 


strained laver beam wayeguide absorber (dashed) at location | and 


without (solid), 
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Figure 40.  Nfodal damping factors vs. frequency of the test plate with a 16” con- 


strained laver beam waveguide absorber at location | and without. 
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Figure 41. The test plate with three waveguide absorbers. 
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Figure 42. The driving point frequency response of the test plate with a 20” viscoe- 
lastic (at location 1) and 16” constrained Jayer (at location 2) beam 


waveguide absorber (dashed) and without (solid). 
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Figure 43. Mfodal damping factors vs. frequency of the test plate with a 20” yiseoe- 


lastic ocation tl) and a 16” constrained layer (location 2) beam wave- 


guide absorber and without. 
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Figure 44. The driving point frequency response of the test plate with a 20” viscoe- 
lastic (at location 1), a 16” constrained layer (at location 2) and a 16° 
viscoelastic (at location 3) beam waveguide absorber (dashed) and 


without (solid). 
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Figure 45.  M{odal damping factors ys. frequency of the test plate with a 20” viscoe- 
lastic (location 1), a 16” constrained laver (location 2) and a 6° viscoe- 


latic beam waveguide absorber and without. 
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Figure 46. Mfagnitude impedances of the test plate at location 1 (solid) and of the 


2()” viscoelastic beam at the driving point (dashed). 
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Magnitude impedances of the test plate at location 2 (solid) and of the 


16’constrained layer beam at the driving point (dashed). 
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Figure 48. Magnitude impedances of the test plate at location 3 (solid) and of the 


16” viscoelastic beam at the driving point (dashed). 
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Vv. CEONCEES Cy 


The purpose of this study was to develop two kinds of high damping beam waye- 
guide absorbers for application to naval structures over a wide frequency range 
(LOOHZ-2000HZ) and to evaluate theoretical prediction schemes using various beam 
theories. Comparison between theoretical predictions and experimental results show that 
Timoshenko beam theory and the six order beam theory can be used in the prediction 
of impedances of viscoelastic beam waveguide absorbers and constrained laver beam 
waveguide absorbers, respectively. 

Application of waveguide absorbers to a test plate structure showed that daniping 
of the test structure can be increased significantly with one waveguide absorber over a 
frequency range. The magnitude and the trequency range of darnping increase depend 
on the inpedance of the waveguide absorber and on the location. The prediction of 
damping increase due to an attached waveguide absorber. equation (2.7). does not show 
the effect of waveguide absorber completely since this equation does not include the 
structure vibration reduction effects due to the waveguide absorber. Since the contrib- 
ution of each waveguide absorber to the damping of the test plate ts additive. the modal 
damping values of all the vibration modes of the test plate ina wide frequency range can 


be increased using multiple waveguide absorbers. 
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VI. RECOMMENDATIONS 


From the results of present study following are recommended for future studies on 
waveguide absorbers: 

1. Investigation of the effect of a waveguide absorber on the reduction in vibration 
of the original structure and derivation of the relation between the impedances and the 
mode shape of the structure, 

2. Theoretical studies on impedance prediction of waveguides of different shape such 
as circular viscoelastic plates and application to damping increase, 

3. Investigation of interaction between waveguide absorber for multiple waveguide 
absorbers applications, and 

4+. Optimization of application of waveguide absorbers in design, number, and 


location. 
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APPENDIX A. 
l. THE IMIPPEDANCE OF THE INFINITE VISCOELASTIC BEANT USTs 


s's eve Je s‘e foves ‘eles evleveslote Seveledeslesles'ecleseVeslese es scevetlsteseveve 
FR FDC ED IDV ED ON FY ED FU GL CB FE FX CHO SHEN HDT FT FH FD SL OH OD GD 


od 


THE THREE REDUCE PROGRAMS WERE USED FOR EVALUATING THE IMPEDANGIIS 
OF THE INFINITE VISCOELASTIC BEAM USING BERNOULLI-EULER BEAM * 
THEORY: 5 


e 
fawawaata ale alan'an's wn wanlenenlawadasatacdaviantsa'svanian’anian wvlanaeclanwasclan seaniavlatawowen’s daclan'en srlan wala we woul ale we ale ween en an's wl =e ate alae ee an en an ‘se enlawen's 
SoS CD FY eB FV GEV. sR ID FER SE CS 6H FR IR CR CD OD SHR ER EDR CR CH 4H SR OD GB Pe EE ER ED FD OD CD SR OVER CD FH CR FRM GR FR FY FY HR FR FB EDV ED FR Gd FB CR SB HR SL FH CR FR CH FH GR ER oH CR FL FH ed CD 


THE FIRST REDUCE PROGRAM DEFINED TRANSVERSE FORCE AND TWO BOUNDARY 
CONDITIONS. 


2a Y 1S EXPCI KI) YS ENP Cla), 


PP?=DEC YS. 

MMS==F De CPP a): 

FF: =-DF(MM,X); 

YO: =SUB(X=0 , YY) -Y; 

PO: =SUB( X=0, PP) -P; 

Bye. 

C 

C THE SECOND REDUCE PROGRAM SOLVED CONSTANT VARIABLES (Y1 AND Y3)3 
C 

YO: =Y 1+Y 3-7: 

POS=F lye: Yo -r: 

SOLVE CES Gyo Ve Ue is) 

Bart. 

C 

C THE THIRD REDUCE PROGRAM EVALUATED THE TRANSVERSE FORCE. 


C 


MYO=i 1 ESE Gieikess ) SEN PCI'ho-%): 
Per =E pray... 


FF: =SUB(X=0,FFF); 

BYE; 

Jeverte s'eve s'e evlevetleseveveveve esevevevevss's eveveveVeJeveres'e ves fowtos ‘eve ss'e evesleveveveves's ves'cs's s'e s': fevestes! eveve s'e veves ‘eslesle s'e sveveVeseve 'eTs's72 
vs THIS FORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF TH= INFINITE vs 
vs VISCOELASTIC BEAM USING BERNOULLI-EULER BEAM THEORY. ¥ 
¥ DEFINED VARIABLES; ve 
¥ E COMPLEX YOUNG'S MODULUS OF THE VISCOELASTIC BEAM % 
vs F : FREQUENCY (HZ) 7 
ve FF : TRANSVERSE FORCE (LBF) ¥ 
vs G : COMPLEX SHEAR MODULUS OF THE VISCOELASTIC BEAN vs 
vs GRAV : GRAVITY (336 IN/SEC 

vs IMP : IMPEDANCE OF THE VISCOELASTIC BEAM *s 
v POISS : POISSON'S RATIO “ 
v ROH : DENSITY OF THE VISCOELASTIC BEAM * 
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hoe ee ENS lONSwer THE waeaN CWIDTH, HEIGHT, LENGTH) 


i : AREA MOMENT OF INERTIA e 
XK1, XK2, XK3, KK4 : WAVE NUMBERS ‘: 


XNETA : ENERGY LOSS FACTOR OF THE VISCOELASTIC BEAM 
v1 : IMAGINARY (0,1) 


RPO ent aloe elo cls clea oalaclactaataats ate uf re tatos oral sclectocts ucts cle cluctoctsa! Weegee aloe cl ea a olan ol nice ala clea sate ale cn alsclnaleclretaclontelec nats elsalsataclocloatontsa! 
BCT C Se IC ro 2e TCC US TC TE DETR TS Fe Ee COCCI Te OTC IC OIC Ie TOT Ree Te TOTTI BE TCI Te TE SOT EUE TE TEI TE ETE TE UE TES ETC T ETE TS TOT EE TONER 


C 


Clara 


ao 2G 


Grea 


COMPLEX V1,G,E,KB1,KB2,XK1,XK2,XK3,XK4,P1,P3,A1,A2,A3,A4,FF, IMP 
REAL F,POISS,GRAV,PI,XB,XH,XG, ROH, XI 
REAL W,XLOG,XNETA,XA,XCA,XSA,KC,IMPR, IMPI 


OPEN(UNIT=15 ,FILE='IFI') 
CONSTANTS 


F=10. 0 
POISS=0. 45 
GRAV=386. 0 

PI=3. 1415926 
V1=CMPLX(0.0,1. 0) 


DIMENSIONS 


XB=1. 0 
XH=0. 35 


PROPERITIES 


ROH=0. 05505/GRAV 
AI=(XB*XH*"*3 ) /12. 0 


FOUR DIFFERENT WAVE NUMBERS 


DO 100 I=1,100 
W=2*PISF 
XLOG=ABS(LOG10(F/62)) 
XNETA=0. 65*EXP( -0. 52732**XLOG**1. 956) 
XNETA=0. 2 
XNETA=0. 5 
XG=0. 00002503*F**3-0, 1752*F%*24+457. 5883*F+29280 
XG=200000 
G=XG**CMPLX( 1. 0,XNETA) 
E=2*G**( 1+POISS) 
XKA=ATAN( XNETA) /4. 0 
KCA=COS( -XA) 
XSA=SIN( -XA) 
KC=( SORT( 1+XNETA*=*2 ) )7"(-0. 25) 
KB1=-CSQRT( ROH*XB*XH*W**2 /(E*XI)) 
KB2= CSORT( ROH*XB*XH*W2 /( E*XI)) 
KK1= CSQRT(KB1)**KC*CMPLX( XCA,XSA) 
XK2=-CSQRT(KB1 )*KC*CMPLX(XCA,XSA) 
XK3=-CSQRT( KB2 )*KC**CMPLX( XGA, XSA) 
XK4= CSQRT(KB2)**KC*CMPLX(XCA,XSA) 
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C DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 


Pi=V1*XK1 

P2=V1*XK2 

P3=V1"*XK3 

PG=V1*XK4 

Al= P3/(P3-P1) 

AZ=0 

AS==F 1 / Ceo -F lo 

A4=0 

FR=-V1IV"E*XIT* (ALS XK 1 34+A3"% G7" 3) 


IMPEDANCES 


aoa © 


IMP=2*FF/(V1*W) 
IMPR=REAL( IMP) 
IMPI=AIMAG( IMP) 
WRITECI > 1000) , THERE aie I 
F=F+20. 0 
100 CONTINUE 
TOGO hORMARCLS 09.5524 ,012 055.2. (pa 
CLOSE(UNIT=15 ) 
SLOP 
END 


2. JHE INfPEDANCE OF THE FINITE VISCOELASTIC BEAT USING 


BERNOULLI-EULER BEAM THEORY 

YY: =Y1"EXP( I°*K1*X)+Y2*EXP( [*K2*X )+Y3"EXPC 1*K3*X)+YS EXEC lho ae 
Pes =DECiyY sx); 

MM: =-ET*DFCPP,X); 

FF: =-DFCMM,X); 

YO: =SUBCX=0, YY) -Y; 

PO? =SUBUA=0. Pe -F: 

ML: =SUBC X=L,MM); 

FL: =SU3(X=L,FF); 

BYE; 

C 

C THE SECOND REDUCE PROGRAM SOLVED CONSTANE VARIABLES “( Y ey Zee 
S AND Y4). 

C 

YOu ey l+yY2rY 374-7. 

PQ: =P iY 1EP 2-243" Y3tra 4 4-P: 

ML: =M1*Yl+M2**Y2Z4+M3**Y3+M4"*Y4; 

FL: =Fi*Y1ItR2*Y2+F3*Y3+R4"Y4, 
SOLVECLSTCYO, POsiL Fi} yey 2 oe eee 

BLS: 

c 

C THE THIRD REDUCE PROGRAM EVALUATED THE TRANSVERSE FORCE. 

C 

YY: =YIVEXP( 1*K1*K)+Y2*EXP( I*K2*K DEY 3*"EXP( E*K 36) eee 
PRP eSB ie (Oty Yorker) 

FF: =SUB( X=0 , FFF); 

BYE: 


Sesesesccccesesesevevesevesescacsdescvevesesesesevcvlevesesledevesievevedcvevedsdeslevevesdesedsdedtevevevevevede te seve levee dete Fe Fete Te eee Be Se 


* 


¥ UiicerORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF THE FINITE 
w ViSCOERASTIC BEAM USING BERNOULLI-EULER BEAM THEORY. 

ae 

¥ DEFINED Veh EeBRES: 

% E : COMPLEX YOUNG'S MODULUS OF THE VISCOELASTIC BEAM 
% F PERE OURNCY (Ha) 

¥ FF PORANS VERSE PORCH CLE) 

* G > COMPLEX SHEAR MOBULUS OF THE VISCOELASTIC aa 

¥ GRAV se ORAVI EY (3a IN7SEC 

% LHP ene EUANCE OF THewy BseCORLASTIC BEAM 

* FOSS : POISSON S RATIO 

¥ ROU POENSTEY*Or THE VISCOELASTIC BEAH 

%* eee ten OL ENS IONS OF THE BEAM (WIDTH, HEIGHT, LENGTH) 
% Xl AREA MOMENT OF INERTIA 

¥ Pian enh, AKG > WAVE NUMBERS 

% XNETA > ENERGY LOSS FACTOR OF THE VISCOELASTIC BEAM 

%* V1 SeMAGINARY CO, 1) 

ae 

Ne ee C10 CoC IE WE Ne ONE eo 207s te aCe COC ICTS IN Te eC IC TE ee CUE TT US IE TR INTC IE ROSS SEC EN GS ICDS TS UR TOTS U OTE Oe CTE TE 
C 


COMPLEX V1,G,E,KB1,KB2,XK1,XK2,%K3,XK4,P1,P2,P3,P4 
COMPLEX M1,M2,M3,M4,F1,F2,F3,F4,DEN,A1,A2,A3,A4,FF,IMP 

REAL F,POISS,GRAV,P1,XB,XL,XH,XG,ROH,XI,W,XLOG,XNETA,XA,XCA,XSA 
REAL KC,KR1,KR2,KR3,KR4,KI1,KI2,KI3,KI4, IMPR, IMPI 


C 
OPEN(UNIT=15 , FILE='FFI' ) 
C 
C CONSTANTS 
C 
F=10.0 
POISS=0. 45 
GRAV=386. 0 
PI=3. 1415926 
VI=CMPLX(0. 0,1. 0) 
c 
C DIMENSIONS 
C 
heal, (0 
=o. 0 
XH=0. 35 
c 
G PROPERITIES 
C 
ROH=0. 05505/GRAV 
C 
C FOUR DIFFERENT WAVE NUMBERS 
C 
XI=(XB*XH*"*3)/12. 0 
DOeL09 I=1, 100 
W=2"*PI“F 
XLOG=ABS( LOG10( F/62)) 
XNETA=0. 65**EXP(-0. 52732**XLOG***1. 956) 
c XNETA=0. 2 
¢ XNETA=0. 5 


US 


Ca Gar Ca 


Rm & 


Q QR © a 


XG=0. 00002503*F**3-0. 1752*F**24457. 5883. ae 
xG=200000 

G=XG*CMPLX(1.0,XNETA) 

b=2 Ghee hss) 
KA=ATANCXNETA) /4. 0 

XCA=COS( -XA) 

XSA=SIN( -XA) 

KO=( SOR TOieaNE TAs*<2)) CC =Oneae) 
KB1=-CSQRT( ROH*XB* XH Were? /CE*XI ) ) 
KB2= CSQRT( ROH*XB*xH*We**2 /CE*XT) ) 
XK1= CSQRTCKB1)*KC*CMNPLXCXCA, XSA) 
XK 2=-CSQRT(KB1)*KC*CMPLXCXCA,XSA) 
XK3=-CSQRT(KB2)*KC*CNPLX( XCA, XSA) 
XKS= CSQRTCKB2)*KC*CMPLACXCA, XSA) 
KRI=REAL(C XK1) 

KR2=REAL( XK2 ) 

KR3=REALC XK3 ) 

KRG=REAL( XK4) 

KI1=AIMAG(XK1) 

KI2=AINMAG( XK2) 

KI3=AIMAG( XK3 ) 

KIG=AIMAG( XkK4 ) 


DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 


P1=V1"*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 

MISEXP( -XL*KI1)*( COSC XL*KR1)+V 1*SIN( AGRI) AR 2 

M2=EXP( -XL*KI2)*( COSC XL*KR2)+V1*SINC XL*KR2) )*XK2** 2 

M3=EXP( -XL*kI3)*( COS( XL*KR3)+V1*S INC XL*KR3 ) )*XK3%*2 

M4=EXP( -XL*KI4)*( COSC XL*KR4.)+V1*S INC XL*KR§S ) ) *XK4**2 

FI=EXP( -XL*KI1)*( COSC XL*KR1I)+V1*SIN( XE=KR1) eR eas 

F2=EXP( -XL*KI2)*( COS( XL*KR2)4+V1°SINCXL*KR2) )#XK24%3 

F3=EXP( -XL*KI3) *C( COSC XL*KR3 )+V1*S INC XL*KR3) )*XK3**3 

F4=EXP( -XL*KI4)*( COSC XL*KR4)+V1*SIN( XL*KRS) )*XK4 7" 3 

DEN=( F4*M3*P2-F4*M3*P1-F4*M2*P3+F4*N2*PI+F4*M1l*P5-F4tisr 2 
“F3*M4P24+F3*M4*P1+F3"M2*P4<-F 37 M2"PI-F3*i1*Po7k oles 
+F2°M4*P3-F2°H4*P1-F2¢N3*P4th 2°U3"0 18 2h Pao b 2 
“FIsN4*P3+F1*M4*P2+F 1*M3*P4-F 1*M3%P2-F 1*M2*P4+F 1*M2*P3) 

Al=( F4*M3*P2-F4"N2"P3-F3"N4°P2+F3*M2*P4-rk 2-M4<P3-F 23) (eo 
/DEN 

A2=(-F4*M3*P1I+F4*M1*P3+F3*N4*P 1 -F3"M1*P4-F lia“ Po+i sie 
/DEN 

A3=( F4*M2*P1-F4*M1*P2-F2*MG*PI+F2*M1“P4+k l“M4“P2-F l“M2"Pa) 
/DEN 

A&=( -F3*M2*PI+F3*M1*P2+F2°M3"P1-F2"MisP3-F 1 tis Ba ee 
/DEN 

FR=-V1ISE*XI*S CALS XK 1S 34+A2"XK254 344A 3° XK 35% 34+AG XRG 3 ) 


IMPEDANCES 
IMP=2*FF/(V1*"*W) 


IMPR=REAL( IMP) 
IMPI=AINMAG( IMP) 


SO 


WRITE( 15,1000) F,IMPR,IMPI 
F=F+20. 0 
100 CONTINUE 
1000 FORMAT(1X,D9. 3,2X,D12.5,2X,D12. 5) 
CLOSE( UNIT=15) 
STOP 
END 


3. THE IMPEDANCE OF THE FINITE VISCOELASTIC BEAN: USING 
eNO ie ae ee 


Tevedeveteveve sede Fete de sete te Fe Fe Fete FETS ETE TE TCT SESE EE TENE SC EVO G ETE TEETER TERETE EE TOE TEETER TET TET TE TE TE TEE BE TE TEETER EE TEESE 
we $3 
* THEY THREE REDUCE PROGRAMS WERE USED FOR EVALUATING THE IMPEDANCE * 
* Of tiie visCORGAS! LC BEAM USING TIIOSHENKO BEAN THEORY. * 
¥ ve 
EIS ISTO TCI CRETE IE IC IE TEI OS ITE IE IETS ICTR UE IETS ICES 
C 

C THE FIRST REDUCE PROGRAM DEFINED SHEAR FORCE AND FOUR BCUNDARY 
C CONDITIONS. 

C 

meee 1 BE XPC LK IX FY ZSEXP( I*K2*X )+Y3"EXPC I1*KR3"X )+Y¥4"°EXPC I*K4X) 5 
Ped: =R1*Y1*EXPCL*K1*X)+R2*°Y2*EXP(C I*K2*X); 

PP2: =R3*Y3*"EXPCI*K3*K)+R4*Y4*"EXP( [*K4*X) 5 

igen —PP1+PP2; 

MM: =DF(PP,X); 

Bo: —~-(DFCYY,X)-PP); 

YO: =SUB( X=0 , YY) -Y; 

eee—SUBC(X=0,PP)-P; 

ML: =SUBC X=L,MM); 

SL: =SUBCX=L,S5S); 

Bie; 

C 

C tite onsOND REDUCE PROGRAM SOLVED CONSTANT VARIABLES (Y1, Y2, 

C Y3 AND Y4). 

C 

mee —YItY2+Y3+Y4-Y; 

Peek 15Y 1 +R2*Y 2t+R3*Y3+R4*Y4 -P; 

Pie HY 1+AN2*Y24+M3°°Y34+M4* 4; 

ee» 1 Y¥14+S2*Y24+83*13+54*Y4; 

SOLVE( LST(Y0,P0,ML,SL) ,Y1,¥2,Y3,Y4); 

PL ES 

C 

C THE THIRD REDUCE PROGRAM EVALUATED THE SHEAR FORCE. 

C 

Y¥: =Y1*EXP( I*K1*X)+Y¥2*EXPC I*K2**X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 
See LY LE XPCT*K1*XK)+R2*Y2*EXP( I*K2*xX); 

PP2: =R3*Y3*EXP( 1*K3*X)+R4°Y4*EXP( I*K4*X); 

eee =F P1+PP2; 

9558: =-KA*GG*AA*(DFCYY ,X)-PP); 

SS: =SUB( X=0,SSS); 

Bye : 


? 
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¥ vc 
w THIS FORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF THE FINITE Ww 
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E COMPLEX YOUNG'S MODULUS OF THE VISCOELASTIC BEAM 
F : FREQUENCY (HZ) 
(5 : COMPLEX SHEAR MODULUS OF THE VisGer bee ee 
GRAV = GRAVITY C32 ba hy see 
IMP ; IMPEDANCE OF THE VISCOELASTIC BEAr 
KAR ; SHAPE “CP CROs Shel TON 
FOTSS : POISSON'S RATIO 
ROH : DENSITY OF THE VISCOELASTIC BEAM 
55 : SHEAR FORCE (LBF) 
AB, XH,XL : DIMENSIONS OF THE BEAM (WIDTH, HEIGH! Eh cae 
XI AREA MOMENT OF INERTIA 
hie XK2, XK3, XK4 : WAVE NUMBERS 
XNETA : ENERGY LOSS FACTOR OF THE VISCOELASII@ eal 
V1 : IMAGINARY (0,1) 
sevededededesedsteveseve ted leds tess evs Vese ee Tess Tea eee se Teese seve Fe eee se ae 08 seas eh ae ee ae Oe ae ae OC 7G Teas on S709 On cea 


C2 2-0) 


C2 €])'C2-0) 


C2 Cc) Ca 


VISCOELASTIC BEAM USING TIMOSHENKO BEAM THEORY. 


DEFINED ee 


COMPLEX V1,G,E,XKN1,XKN2,XKD,XK1,XK2,XK3,“K4,R1,R2,R3,R4 
COMPLEX NM1,M2,M3,M4,$1,$2,$3,S4,DEN,A1,A2,A3,A4,SS, IMP 
REAL F,KAR,POISS,GRAV,PI,XB,XL,¥H,G,ROH,XI,W,XLOG, XNETA 
REAL KR1,KR2,KR3,KR4,KI1,K12,KI3,KI4, IMPR, IMPI 


OPEN( UNIT=15 , FILE='TIM' ) 
CONSTANTS 


F=10. 0 

KAR=0. 83 
POTSS=0. 45 
GRAV=386. 0 

PI=3. 1415926 
VI=CMPLX( 0. 0,1. 0) 


DIMENSIONS 


XB=1. 0 
XL=8. 0 
AH=0. 35 


PROPERITIES 


ROH=0. 05505/GRAV 
XI=(CXB*XH**3)/12.0 


FOUR DIFFERENT WAVE NUMBERS 


DO 100 I=1,100 
W=2"PI“F 
XLOG=ABS(LOG1O(F/62)) 
XNETA=0. 65*EXP(-0. 52732*XLOG**1. 956) 
XNETA=0. 2 
XNETA=0. 5 


— Q © ee 


KRG=O0300002503*F **3-0. 1752*F* 24457. 5863*F+29280 

XG=200000 

G=XG*CMPLX(1.0,XNETA) 

E=2*G*(1+POISS ) 

XKN 1=W****2**( KAR*G+E ) /ROH 

XKN2= SQRTCW**4*( (KAR*G+E ) /ROH) **2 
-4*KAR*G*E*( WG -KAR*G*XB*XKH*Ws* 2 /( ROH* XI) ) 
/CROH***2 ) ) 

XKD=2**KAR*G*E /( ROH***2 ) 

XK1= CSQRT( (XKN1+XKN2) /XKD) 

XK2=-CSQRT((XKN1+XKN2) /XKD ) 

XK3= CSQRT(( XKN1-XKN2) /XKD) 

XK4=-CSOQRT( ( XKN1-XKN2) /XKD) 

KRI=REAL( XK1) 

KR2=REAL(C XK2) 

KR3=REALC XK3 ) 

KR4=REAL( XK4 ) 

KI1=AIMAG( XK1) 

KI2=AIMAG( XK2) 

KI3=AIMAG( XK3) 

KI4=AIMAG( XK4) 


DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 


R1=( ROH*Ws** 2 -KAR*G*XK 12 ) /( V1" KAR*G*XK1) 
2=( ROH* We" 2 -KAR*G*XK2**2 ) /( V1KAR*G*XK2 ) 
R3=( ROH*W2 2 - KAR*G*XK3%2 ) /( V1*KAR*G*XK3 ) 
R4=( ROH*W*7*2 -KAR*GXK4* 2) /( V1KAR*G*XK4 ) 
M1I=EXP( -XL**KI1)*(COS(XL*KR1)+V1*SIN( XL*KR1) )*V1"°XK1*R1 
M2=EXP( -XL*KI2)*(COS(XL*KR2)+V 1*SIN( XL*KR2 ) )*V1*XK2*R2 
M3=EXP( -XL*KI3)**( COS(XL*KR3) +V1*SIN( XL*KR3) )*V1*XK3*R3 
M4=EXP( -XL*KI4)*( COS( XL*KR4)+V1*SIN( XL*KR4 ) )*V1*XK4*R4 
S1=EXP( -XL*KI1)*(COS(XL*KR1)+V1*SIN(XL*KR1) )*(R1-V1"XK1) 
S2=EXP( -XL*K1I2)*( COS( XL*KR2)+V1*SIN(XL*KR2) )**(R2-V1*XK2) 
S3=EXP( -XL*KI3)*(COS( XL*KR3)+V1*S IN( XL*KR3) )*(R3-V1*XK3) 
S4=EXP( -XL*KI4)**(COS( XL*KR4)+V1*SIN( XL*KR4) )*( R4-V1*XK4 ) 
DEN=( $47*M3**R2 -S4**M3**R1-S4*M2*R34+S42M2*R1+S47M 1°R3-S47M1*R2 
= S3*M4*R2+5 3°M4GR 145 3%M2"*R4-S32M2"R 1 -S3*M 1 R4+S 32M LR2 
+52*M4"*R3-S2*M4*R1 -S2*M3*R44+S52*M3*R1+52*M1*R4-S2*M1*R3 
=S1M4°*R34+S 1°M4*R2+5 1M 3**RG -S 1M 3*R2-S 12M2*R4+S 1M2*R3) 
Al=( S4**M3**R2-S4%M2*R3-S3%M4*R2+S 3*M2*R4+S 2*M4*R3 -S2*M3*R4) 
/DEN 
A2=( -S4**M3*R1+S4°M1*R3+S3%*M4*R 1 -S3**M1*R4-S1*M4G*R3+5 1°M3*R4) 
/DEN 
A3=( S4*M2*R1-S4*M1*R2-S27M4*R1+S 2M 1R4+S 19M4G*R2 -S 1*M2*R4) 
/DEN 
A4=( -S3%M2*R14+S3*M1*R2I+S 27M 3*R 1 -S 2M 1*R3 -S1*M3*R24S 1°M2*R3) 
/DEN 
SS=KAR*G*XB*XH*( Alt*(R1-V1*XK1)+A2%(R2-V1*XK2 )+A37*( R3-V1*XK3) 
+A4%(R4-V1*XK4) ) 


IMPEDANCES 


IMP=2*SS/(CV1*W) 
IMPR=REAL( IMP) 
IHMPI=AIMAG( IMP) 
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WRITE( 15 , 1000 ERS Mee 
F=F-+20. 0 
100 CONTINUE 
1000 FORMAT(C 1X, Don o62 Dio 2, ore 
CLOSECUNIT=15 ) 
[LOE 
END 


4. THE INIPEDANCE OF THE FINITE CONSTRAINED LAYER BEAN] 
eae Pee Se Oe Cae Ts 


‘ess sles‘ aoe ase eaee £8 2% 2 eat ma Yess se ae ‘ole s< a2 ae oes ales Je on al. aye Cs we aa olen al. Py as ata “a ae ate oa ae aus 
sesicss Sesesses eae ee eae FB se ee Fd FR ED ERS CR FX ER FE ER FT CH FYT FE FX EB ENGR ER GE eS FB ED ER 


¥ THE THREE REDUCE PROGRANS WERE USED FOR EVALUATING THE [MPEDANGEs 

* OF THE INFINITE VISCOELASTIC BEAM USING BERNOULLI-EULER BEAM “ 
THEORY. vs 

3 7% 

Yetovesedevedevededevevede vedevedese ds Kesetete se Ke seteveveseveve seve se see se Begs Wee Ue Ieee IC ee IS DCIS DS ae ae ee 0s oe IS ICTS ye oe OCS oe oe 

C 

C THE FIRST REDUCE PROGRAM DEFINED TRANSVERSE FORCE AND TWO RCUMDARY 

G CONDITIONS. 

G 

YY1i; =Y1"EXPC I1*K1*X) +Y2*EXP( I*K2*X); 

Y¥2: =¥3*EXPCI*KS> kya E PCI Kak). 

YY3? =Y5"EAPCIKS 6 Be ho Or 

UU12 SRY ISEB K PCTS iees )-R Zao 2c BaP aco Zeck) 

UU2: =R3*Y3"EXPC I*K3*X)+R4*Y4*“EXPCI*K4*X); 

UUS; =Ro“ YS" EXPE KS*X) ERO YOr ERE take). 

YY Ye SV evry oe 

UUU: =UU14+UU2+UU3; 

YO: =SUBC X=0, YYY)-Y; 

PO: =SUBCX=0 , DFC YYY,X) )-P; 

U0: =SUB( X=0 , UUU) -U 

SL: =SUB( X=L, ( PP*DF(YYY, X,3)+QQ*DF( YYY, X)+RR*UUU)); 

ML: =SUBC X=, CPP*DF(YYY,A,2)roo* OF CUUUL ee 

NL: =SUB(X=L, (TT*DFCUUU,X))); 


Bes 

C 

C THE SECOND REDUCE PROGRAM SOLVED CONSTANT VARIABLES (Y4, Y¥5, 
C AND Y6). 

C 


VOP = 71 ey oer ey ooo 

POVS=P ie yo (OEP or oe ee ot oo. 
UO? =R i-s¥ [4R2Y2FRS4Y STR 1S thou: 
SL: =Si"Y1i+S2* 253" arse Yor55"Ystoar yo: 
ML: =M1I*Y14M2*Y2+H3*Y34N4“Y4+H5*Y54+M6"1 6; 
NOU: =N1*¥ 1-4N2*°Y 24N3*YS4+N4°Y44N5“Y5-+N6* YO; 
SOLVECEST( OLE iy 4 uo. 

BCE: 


THE THIRD REDUCE PROGRAM SOLVED VARIABLES (¥1, Y2 AND Y3) USING 
Y4, YS AND Y6. THE THIRD REDUCE PROGRAM EVALUATED THE SHEAR PORGEs 


Co, GasCa, Ca 


V4: =B1'¥14+C1*Y2+D1*Y3; 
Y5: =B2**Y1+C2"*Y2+D2*Y3: 
Y6: =B3*Y14C3*Y2+D3"Y3; 
YO: =Y1+Y24+Y34Y4+Y5+Y6-Y; 


PO: =P1*Y1+P2*Y¥2+P3*Y3+P4*Y4+P5*Y5+P6*Y6-P; 
UO: =RI*V14+R2*Y2+R3*Y3+R47° Y44R5*Y5+RO"Y6-U; 
SOLVE(LST(Y0,P0,U0),Y1,Y2,Y3); 

YYY 1: =Y1*EXP( [°K 1*X)+Y 2" EXP(1*K2*X); 

YYY2: =Y3*EXP( [*K3*X)+Y4*EXP( I*K4*X); 

YYY3: =YS*EXP( [*K5**X)+Y6%EXP( [*K6"X); 

YYYY: =YYY1+YYY2+YYY3; 

UUUL1: =R1* YL EXP( TK 1X )4+R2"Y 2" EXP(I*K2*X); 
UUU2: =R3*Y3*EXP( 1#K3*X)+R4*Y4"EXP( T*K4#X),; 
UUU3: =R5*Y5*EXP( I*K5*X) +R6*Y6*EXP( [*K6"X); 
UUUU: =UUU1+UUU2+UUU3; 

SSS: =PP*DF( YYYY,X, 3) tOQ*DF(YYYY , X)+RR*UUUU; 
SS: =SUB(X=0,SSS); 


DE, 
Sevesicsevesesesedesclesedevesevevesedevedevevevevevevedevevevededeteve deve te deve ds ete ve Fe Fe dese cece te dete TS Ge LOGE ET CEE ETC STE TS TEE GE TE 
ve 
7% THIS FORTRAN PROGRAM DEVELOPED THE IMPEDANCE OF THE CONSTRAINED » ¥ 
¥ AYER BEAM USING THE 6TH ORDEK BEAM THEORY. 
Jc 
¥ DEFINED VARIABLES; 
¥ El, E3 : COMPLEX YOUNG'S MODULUS OF THE CONSTRAINED LAYER 
16 F : FREQUENCY (HZ) * 
¥ G PVCOME LE SHEAR MGDULUS OF MHiEavesCOELASTIC * 
* GRAV SeGRAV EY (386 IN/SEC 72) * 
¥ elo weeeae  Ght OF THE CONSTRAINED LAYER * 
¥ H2 SHE LoOlie Ors tHeey ISsCOEnaslic vw 
¥ ie 12 os AREA MOMENT OF INERTIA if 
* rte Sere DANCENOue THE] CONSTRAINED LAYER BRA * 
* KAR Ss SHAPE OP SCRGS5 SECTION * 
* POISS : POISSON'S RATIO 
% ROH I SSUENSI DY OF THE GONSTRAINED LAYER 
% ROH2 PeUbho bY SOE THE VISCOELASTIC 
¥ oo Pe onn AR EFORCE, (EBF ) 
# hp, Au =| DIMENSIONS OF THE BEAM (WIDTH, LENGTH) 
* XK1, XK2, XK3, XK4, XK5, XK6 : WAVE NUMBERS 
ve XNETA metwERoveL0o> FACTOR OF THE VISCOELASTIC BEAM - 
¥ V1 : IMAGINARY (0,1) * 
ve 3 
Terese devedese dees He seve vowels seve we ess ae Tete eve ve Tes Toe Te Ne Te Te Tete PETE TS Se ETE TS ES ETS NE TE TELE TOTS TS OTE TS ETE PETE TE ESTEE TS 
C 
COMPLEX*32 V1,G,BB,P1,P2,P3,P4,P5,P6,R1,R2,R3,R4,R5,R6,Q0,RR 
COMPLEX**32 $1,S2,53,S4,55,S6,M1,M2,M3,N4,M5,M6,N1,N2,N3,N4,N5 ,N6 
COMPLEX*32 B1,C1,D1,B2,C2,D2,B3,C3,D3,DEN1,DEN2,A1,A2,A43 
COMPLEX B,D,XK1,XK2 ,.XK3,XK4,XK5 ,¥K6 ,XXK1,XXK2,XXK3, IMP 
emt wowace GRAVY XB. AL,XG,H1,H2,H3,E1,53,11,13,DD,ROH1,ROH2 
REAL*16 CC,XM,T1(30) ,T2(30) ,T3(30) ,01(30) ,Q2(30) ,Q3(30) 
Weegee AbOG W,2,C, BK, BI.CR,DR,DL,FT1,FOL,DITI,DITO1.DOO1 ,DQTL 
REAL*16 FT2,FQ2,DTT2 ,DTQ2 ,DQQ2 , DOT2,FT3,FQ3,DTT3 ,DIQ3 , DQQ3 ,DOT3 
Revue to MRRIKKR2,XKR3,XK11,XK12 ,XKI3,KR1,KR2,KR3,KR4,KR5 ,KR6 
REAL*16 KI1,KI2,KI3,KI4,KI5,KI6,PP,S3,TT,IMPR,IMPI 
REAL ANETA 
C 
OPEN( UNIT=15, FILE='SAN' ) 
C 


C2. 


aqna 


<2 2 02 


Ci Ga Ga Co 


CONSTANTS 


F=2010 

PI=3. 1415926 
GRAV=386. 0 
V1=CMPLX(0. 0,1. 0) 


DIMENSIONS 


XB=1. 0 
XL=8. 0 
H1=0. 0625 
H2=0. 065 
H3=0. 0625 


PROFERITIES 


E1=10000000. 0 
E3=10000000. 0 

L1=XBeH L"3 /12. 
13=XB"H3%3 /12. 

DD=E 1*11+E3*13 

ROH1=0. 09832/GRAV 

ROH2=0. 03663/GRAV 

Cea) 2oeeae 

XM=XB**( (H1+H3 )**ROH1+H2*ROH2) 


INITIAL GUESS WAVE NUMBERS 


T1(1)=-4. 31911862 
T2(1)=-0. 05139608 
T3(1)= 3. 77697017 
Q1(1)=-0. 33835599 
Q2(1)=-0. 05104052 
Q3(1)=-0. 27669665 


FIND SIX DIFFERENT WAVE NUMBERS USING NEWTON'S METHOD 
SIX DIFFERENT WAVE NUMBERS ARE XK1,XK2,XK3,XK4,XK5 ,XK6 


DO 100 I=1,100 
W=2"PI*F 
XLOG=ABS( LOG10(F/800) ) 
XNETA=1. 1*EXP( -0. 3906**XLOG**1. 53) 
XNETA=0. 5 
XNETA=1L. 0 
XG=-0. 00001258*F**2+0, 2472"*F+74, 988 
XG=1000 
G=KG*:CMPLX( 1. 0,XNETA) 
BB=G"( E1*H1+E3*"H3) /(H2*E1*H1*E3*H3) 
Z=CC2*E 1S ULE 3°°H3 / (DD ( E1*H1+-E 37913 )) 
B=BB*(1. 0+2) 
C=-XM*W*2 /DD 
D=~-BB*XM*We"2 /DD 
BR=REAL(B) 
BI=AIMAG(B) 
CR=C 


200 


DR=REAL(D) 
DI=AIMAG(D) 
DOMZ00" J=1 51> 
Pl Eedie ore) ieee = oO J)e2-2*B TOL J)+tCR TiC J) 
+( -BR*Q1( J)**2+DR ) 
FQ1=Q1( J) **3+BI%Q1( J) 24( -3*T1( J) 2 -2*BR*T1(J)-CR)*Q1(J) 
rca teimlae 21). } 
Dimes hi J) 22 BR TC) ) 
+( -3%Q1( J) 2 -2**BI*Q1(J)+CR) 
PO eo oes mR econ) )-2*BIisT1( J) 
DQQ 1=3°*Q li J ) Ws PEIKB T*Q 1( J ) 
eo idee = 2 BRT 1 iio) 
DQT1=( -6*QO1( J) -2*BI)*T1C J) -2*BR*O1( J) 
ed =e ee Ceo m  piOl- Fl iy CDT Vi DO01 -bi01- DOT 1) 
Oe On) Older tiie l-POly / (Oli l-DpOOL-pio0l* DoT 1 ) 
Pi ee oR) 2+ oe OZ( J j= 2B i212 J) 
+( -BR*Q2( J)%*2+DR) 
BOZO Je opt O7 (0) 2 = 37 2( J 2-2 BRT 2( J) CR )*Q2( J) 
+ =D TG) 2 =) 
Dig 2=3°12( J es 22 Beet) 
Poe) Ady 2 = 7°°R 1-02 (J )FCR ) 
WMO Z—C= Or Ca Zine Oe Cie 20° B 2 (J) 
DQQ2=3*Q2( J) **2+2*BI*Q2( J) 
+( -3°%°T2( J) "2 -2*BR*T2(J)-CR) 
DOT2=( ~6*02( J) =2"Bi)*22¢ J) -2*BR*Q2( J) 
DAC =F) = NO02Z eT 2=DT02* FOZ) /_ DT12*=DQ02-DI02*DOT2 ) 
O7Glt 1) 0203 = ( -WOtZ Fb PZT 202 )/{ DTT 20002 -DI02*DOT2 ) 
FI3=T3( J)**3+BR*T3( J)**2+( -3"03( J)**2-2*BI*Q3(J)+CR)*T3(J) 
+( ~BR*Q3( J )**2+DR) 
FQ3=03( J)**3+BI*Q3( J)**2+( -3*T3( J) #*2-2*BR*T3( J) -CR)*Q3( J) 
TAs eee Sa)! e826 2 = (I ) 
DTT3=3*T3( J) **2+2*BR*T3(J) 
C= 30 3( J) 2—-2*B1*03( J )+CR) 
DTQ3=( -6*T3( J) -2*°BR)*Q3( J) -2*BI*T3( J) 
DQQ3=3*Q3( J )**24+2*BI*03( J) 
Ge or SC Jee bed o( J) -CR) 
DQT3=( -6*Q3( J) -2*BI)*T3( J) -2*BR*Q3( J) 
T3(J+1)=T3(J)-€ DQQ3*FT3-DTQ3*FQ3) / (DTT3*DQQ3 -DTQ3*DQT3) 
Q3( J+1)=Q3( J) -( -DQT3*FT3+DTT3*FQ3 ) /(DTT3*DQQ3-DIG3*DQT3) 
CONTINUE 
XKR1= T1(J) 
eh 2—= 2) 
XKR3= T3(J) 
i T= ond) 
XKI2= Q2(J) 
XKI3= Q3(J) 
AAKI=CHPLX(CXKR1, XKI1) 
XXK2=CNPLXCXKR2 , XKI2) 
AXAK3=CMPLX(XKR3 , XKI3) 
XK1= CSQRT(XXK1) 
XK2=-CSQRT( XXK1) 
XK3= CSQRTCXKK2) 
NK4=-CSQRT(XXKK2 ) 
XK5= CSQRTCXXK3) 
XK6=-CSQRT( XXK3 ) 
KRi= REALCXK1) 
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CanGa 1G) 


& 
& 
& 
& 


& 
& 
& 


& 


& 
& 
& 
& 


KR2= REAL( XK2) 
KR3= REAL(XK3) 
KR4= REAL( XK4) 
KRS= REAL( XK3) 
KR6= REAL( XK6) 
KLI=AIMAG( XK1) 
KI 2=AIMAG( XK2) 
KI3=A IMAG( XK3) 
KIG=AIMAG( XK4) 
KIS=AIMAG( XKS) 
KI6=AIMAG( XK6) 


DEFINED CONSTANTS FROM THE THREE REDUCE PROGRAMS 


P1=V1*XKk1 
P2=V1*XK2 
P3=V1*KK3 
P4=V1*XK4 
P5=V1*XK5 
P6=V1*XK6 
R1=V1*BB*Z*DD*XK1/( CXK1**2+BB)*E3**H3*CC ) 
R2=V1*BB“Z*DD*XK2/( (XK2***2+BB)*E3*H3*CC ) 
R3=V1*BB*Z*DD*XK3 /( (XK3**2+BB )*E3*H3*CC ) 
R4=V 1*BB*Z*DD*XK4/( (XK47**2+BB) *E3*H3*CC ) 
R5=V1"*BB*¥Z*DD*XK5/( (XK5*"*2+BB)*E3*H3*CC ) 
R6=V1*BB*Z*DD*KK6/( (XKO7"*2+BB) *E3*H3*CC ) 
PP=DD 
QQ=-DD*BB*Z 
RR=DD*BB*Z*( E1*H1+E3*H3 ) /(CC*E1*H1 ) 
SS=-CC*E3*H3 
TT=E3*H3 
o1l=(-Vie=PP=xR1**34+V1"QQXK IRR RI ee ce) 
*( COSC AL KRID+VI<SiNG tie hha 
5S 2=( -V1*PP*XK2*7*3+V 1*QQ*KK24+RR*R2 ) SEXP( -XL*KI2) 
*( COSC XL*KR2 )+V1*SINC XL*KR2) ) 
S3=( -V1*PP*¥XK3*"3+V 1*QQ*XK3+RR*R3 ) *EXP( -KL*KI3) 
*( COSC XL*KR3 )+V1"*SINCXL*KR3) ) 
S4=( -V1*PP¥KK4**°3+V 1 QQ" XKR4+RR*R4) *EXP( -XL*KI4) 
*( COSC XL*KR4)+V1*"SINCXL*KR4) ) 
S5=( -V1*PP*XK5*°34+V1"QQ*XK5+RR*RS5 ) *EXP( -XL*KIS) 
*#( COSC XL*KR5)+V1"SINCXL*KRS) ) 
S6=( -V1*PP*XK6"*3+V1*QQ*X\K6+RR*R6 )*EXPC( -XL*KI6) 
*( COSC XL*KR6)+V1*°S INC XL*KRG) ) 
MI=(V1*SS2R1"xKL-PP XK 2) Exe G- sich 
*( COS( XL*KR1)+V1*SIN( XL*KR1)) 
M2=(V1*SS"R2=5K2 PP aa 2 2 eb ea») 
*(GOSG.LeER2 TV aol h ina 
M3=(V1*SS*R3*XKK3-PP*XKK3**2 )*EXP(-XL*KI3) 
*( COSC XL*KR3)+V1L*SIN(C XL*KR3) ) 
M4=(V1"SS*R4*KK4-PP*KK4"*2 ) "EXP(C -XL*KIG) 
**( COS( XL*KR4)+V1*SIN(XL*KR4)) 
M5=(V1*SS*R5"XK5 -PP*KK5 "2 ) “EXP( -XL*KI5) 
*( COS( XL*KR5 )+V 1S INC XL*KR5) ) 
M6=(V1*SS*R6*XK6 -PP*XK6"* 2 )*EXP( -XL*KI6) 
*( COSC KL*KR6 )+V1*SINCXL*KR6) ) 
N1=S=V1*TT*RIsXKI*EXP( -XL*KI1) 
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&  %*(COS(XL*KR1)+V1*SINCXL*KR1) ) 

N2=V1*TT#R2"*XK2*EXP( -XL*KI2) 

&  *(COS(XL*KR2)+V1*SIN(XL*KR2) ) 

N3=V 1" TT" RS" XK3*EXP( -XL*KI3) 

&  *(COS(XL*KR3)+V1*SIN( XL*KR3) ) 

NG=VL*TT*R4*XK4G°EXP( -XL*KIG) 

&  %*(COS(XL*KR4)+V1*SIN( XL*KR4)) 

NS=V1*TT*R5*XK5**EXP( -XL*KIS5) 

&  %*(COS(XL*KRS )+V1*SIN( XL*KRS ) ) 

N6=V 1" TT*R6*XK6*ENP( -XL*KI6) 

&  *(COS(XL*KR6)+V1*S IN( XL*KR6) ) 

B1=( -N6*M5%S1+N6"M1*°S5+N5*M6*S 1-NS5*1H1"S6 -N1*M6"S54+N1°N5 °S6 ) 
& /(N6*MS5*S4-N6%*M47*S5 -N5*M6*S44N5°N4*S6TNG*M6%S 5 -NGHM5°S6 ) 
C1=( -N6*M5**S2+N6%*M2*S5 ENS *M67S 2 -M5%*M2"*S6 -N2*M6"S5+N2*M57*S6 ) 
& /(N6*MS*S4-N6*M4*S5 -N5*M6%S44N5*MGS6+NG*M6%S5 -NG*M5 56 ) 
D1=( -N6*N5*S 34+N6*M3*S 5+N5°*M6%S 3 -N5%M3°S6 -N3"N6*S5+N3°*M5*S6 ) 
& /(N6*M5*S4-N6%*M4*S5 -N5*M62°S44N5 MGS 6FN4G*M6%S5 -NGH5 7S6 ) 
B2=(N6%*M47°S 1 -N6*M1"*S4 -NG*M62"S 1L4NG*M 1°S6+4N1°M6%S4 -N1M4*S6 ) 
& /(N6*M5*S4-N6*M4*S5 -N5*M6*S44N5 MGS 64N4*M6%S5 -NG*M5 #56) 
C2=( N6*M42S 2-N6*M2*°S4-NG*M6"S 24NGM2 *S6+N2*M6%S4 -N2*MGIS6 ) 
& /(N6*M5*S4-N6*M4%*S5 -N5*M6*SG4N54MG*S6+NG*M6S5 -NGUMS*S6) 
D2=(N6%M4**S 3 -N6*M37*S4 -NG*M6%*5 34N4%M 32S 6+N3*M6%S4 -N3%M4S6 ) 
& /(N6*M5*S4-N6*M4*S 5 -NS*M6*SG+N5%M4*S 6+NG*M6%S5 -NGUM5*S6 ) 
B3=( -N5%*M42*S 14N5%*M1*SG4NG*M5*S 1 -NG*M1*S5 -N1MS*SGtN1°M4*S5 ) 
& /(N6*M5*S4-N6*M4*S5 -N5*M6*S44N5°M4*S6TNG*M6%S5 -NG*M5 #56 ) 
C3=( -N5*M4*S 2+N5%M27*°S44NG*M5*S 2 -N4°*N2°S5 -N2°M5 *S44N2*M47S5 ) 
& /C(N6*M5*S4-N6*M4*S5 -N5*M6*SGENS “HAS 6FNG*M6%S5 -N4EM527S6 ) 
D3=( -NS*M4*S3+N5*M3%S44+NG%M5 7S 3 -N4G*H3%S5 -N3%M5**SG4N3%M4S 5 ) 
& /CN6*M5*S4=N6%M4*S5 -N5*MG6*S4ENS *M4*S6FNG*M67S5 -N4M5 2°86 ) 
DEN1=( R6*P5 -R5**P6 )*( D3*C2*B 1+D3*C2-D3*B2*C1-D3*B2-C3*D2*B1 

. -C3*D2+C3*B2*D1+C3*B2+B3*D2*C1+B3*D2 

& -B3*C2"D1-B3*C2) 

&  +(R4*P6-R6*P4)*(D3*C2*B1-D3"*B2*C1-D3"*C1+D3*B1 -C3*D2*B1 

e +C3*B2*D1+C3"D1-C3*B1+B3*D2*C1-B3*C2"D1 

& -B3**D14+B3*C1) 

& +(R5*P4-R4*P5 )*(D3*C2**B1-D3*B2*C1-C3*D2"B1+C3*B2*D1 

fe +B3*D2*C1-B3*C2*D1+D2*C1-D2*B1-C2*D1 

& +C2*°B1+B2"D1-B2*C1) 
DEN2=(R3**P6-R6*P3)*( C3*B2+C3**B1+C3-B3*C2-B3*C1-B3) 

&  +(R6*P2-R2*P6)*(D3*B2+D3*B1+D3-B3*D2-B3*D1-B3) 

& +(R1*P6-R6*P1)*(D3*C24+D3"C1+D3-C3"D2-C3*D1-C3) 

&  +(R5*P3-R3*PS )*(C3*B2-B3%C2-C2*B1-C2+B2"C1+B2) 

& +(R2*P5-R5*P2 )*(D3*B2-B3*D2-D2"B1-D2+B2*D1+B2) 

&  +(R5*P1-R1*°PS5)*(D3*C2-C3*D2-D2*C1-D24+C2"*D1+C2) 

& +(R4*P3-R3*P4)*(C3%*B1-B3**C1+C2*B1-B2*C1-C1+B1) 

&  +(R2**P4-R4*P2)**(D3*B1-B3**D1+D2"B1-B2"*D1-D1+B1) 

&  +(R4*P1-R1*P4)*(D3*C1-C3*D1+D2*C1-C2"D1-D1+C1) 

& +(R3*P2-R2*P3)*(B3+B2+B1+1) 

&  +(R1*P3-R3*P1)*(C3+C2+C141) 

& +(R2*P1-R1*P2)**(D3+D2+D141) 

A1l=( R6*( P5*D3*C2-P5*C3*D2+P4*D3*C 1 -P4*C3"*D1-P3*C3+P2*D3) 
& + R5*( -P6*D3*C2+P6*C3*D2+P4*D2*C 1-P4"*C2"D1-P3*C2+P2*D2) 
& + R4*(-P6*D3*C1+P6*C3*D1-P5*D2*C1+P5*C2*D1 -P3*C1+P2*D1) 
& + R3*( P6*C3+P5*C2+P4*C1+P2) 

& + R2*(-P6*D3-P5*D2-P4*D1-P3) )/(DEN1+DEN2) 
A2=( R6*( -P5*D3*B2+P5*B3"D2-P4*D3"B1+P4*B3*D1+P3*B3-P1"D3) 
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ca 


100 
1000 


+ R5*( P6*D3*B2-P6"*B3*D2-PS*D2*B1l+P4*B2“D i obo een 
+ R4*( P6*D3*B1-P6*B3*D1+P5*D2"B1=P5"B2-) eo o ee eee 
+ RO C6" Bo -P> Boar eee 

+ R1I*¥( P6é*D3+P5*D2+P4*D1+tP3S) ya ve hieaeiee 

A3=( R6*( P5*C3**B2-P5*B3*C2+P4*C 3B 1 Pa Bo he ee eee 

+ R5*( -P6*C3*B2+P6*B3*C2+P4"C2*Bl=-PssB2-C l=F 2 bea ee 
+ R4*( -P6*C3*B1+P6*B3*C1-P5*C2*8 ites" BZ Ge 2 Pee ee) 
+ R2°( (Po "Bore se Be+ Pood) 

+ Ri*( =P6"C3=P5*C2-P4"C1-P2) 9 7 (Orn ieee ee 


OH KI > Qo —& 


IMPEDANCES 


IMP=2*( (CR4*RR+V1*QQ*XK4-V1*PPYKK4 73 (D1 A34+C1*A2+B 141) 
+ (RS*RREV1*QQ*XKS -VIsPP*XKS*7°3 (D2 °A 34024 ee 
+ (R6*RR+V1*QQ*XK6-V1*PP*XKG67"3 (D3 A3+C3*A2+B377A1 ) 
+ RR*CR3*A3+R2*A2+R1*A1 ) 
+ V1*QO*(XK3*°AS+KK2*A2+KK1*A1) 
-Vi*PP*( K 3% 3% A SEKI 3 AILK 13 %*A1) y/CV1*W) 

IMPR=REAL( IMP) 

IMPI=AIMAG( IMP) 

WRITEC 15 OCC) PF Turk urd 

F=F-20 

TiC PHT 1 

T2023 

T3G) =15( 3) 

OUGY) or) 

Q2(1)=Q2(J) 

Q3(1)=03(J) 

CONTINUE 

FORMAT( 1X,D9. 3,2X,D12.5,2X,D12. 5) 

CLOSE( UNIT=15) 

STOP 

END 


i? Ce 
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Figure 49. Real part of the 20” viscoelastic beam waveguide absorber impedances 


at the center of the beam. 
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Figure 50. Imaginary part of the 20” viscoelastic beam waveguide absorber impe- 


dances at the center of the beam. 
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Figure 31. The driving point frequency response of the test plate with a !6” yiscoe- 
lastic beam waveguide absorber (dashed) at location | and without 


(solid). 
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fodal damping factors vs. frequency of the test plate with a 106” viseoe- 


lastic beam waveguide absorber at location { and without. 
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The driving point frequency response of the test plate with a 20” con- 
Strained layer beam waveguide absorber (dashed) at location | and 


without (solid). 
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Figure 54.  M{lodal damping factors vs. frequency of the test plate with a 20” con- 


strained layer beam waveguide absorber at location 1 and without. 
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Figure 55. The driving point frequency response of the test plate with a 16” viscoe- 
‘ 


lastic beam waveguide absorber (dashed) at location 2 and without 


(solid). 
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Figure 36.  Mfodal damping factors vs. frequency of the test plate with a 16” viscoe- 


lastic beam waveguide absorber at location 2 and without. 
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Figure 57. The driving point frequency response of the test plate with 2 16” con- 
S SP P p 
strained layer beam waveguide absorber (dashed) at location 2 and 


without (solid). 
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Figure 38. Modal damping factors vs. frequency of the test plate with a 16” cun- 


strained laver beam waveguide absorber at location 2 and without. 
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Figure 59. The driving point frequency response of the test plate with a 20” con- 
strained layer beam waveguide absorber (dashed) at location 2 and 


without (solid). 
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Figure 60.  Nlodal damping factors vs. frequency of the test plate with a 20” con- 


Strained layer beatin waveguide absorber at location 2 and without. 
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The driving point frequency response of the test plate with a 16” viscge- 


lastic beam waveguide absorber (dashed) at location 5 


(solid). 
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Figure 62. Modal damping factors vs. frequency of the test plate with a 16” viscoe- 


lastic beam waveguide absorber at location 3 and without. 
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Figure 63. The driving point frequency response of the test plate with a 16” con- 


Strained layer beam waveguide absorber (dashed) at location 3 and 


without (solid). 
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Figure 64.  fodal dumping factors ys. frequency of the test plate with a 16” con- 


strained laver bearn waveguide absorber at location 3 and without. 
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Figure 65. The driving point frequency response of the test plate with a 20” con- 


~ 


strained layer beam waveguide absorber (dashed) at location 3 and 


without (solid). 
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Figure 66.  Mfodal damping factors ys. frequency of the test plate with a 20” con- 


Strained layer beam waveguide absorber at location 3 and without. 
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